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The possibility to use perturbation theory to systematically improve calculations on circular 
quantum dots is investigated. A few different starting points, including Hartree-Fock, are tested 
and the importance of correlation is discussed. Quantum dots with up to twelve electrons are treated 
and the effects of an external magnetic field are examined. The sums over excited states are carried 
out with a complete finite radial basis set obtained through the use of B-splines. The calculated 
addition energy spectra are compared with experiments and the implications for the filling sequence 
of the third shell are discussed in detail. 
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I. INTRODUCTION 

During the last decade a new field on the border be- 
tween condensed matter physics and atomic physics has 
emerged. Modern semi-conductor techniques allow fab- 
rication of electron quantum confinement devices, called 
quantum dots, containing only a small and controllable 
number of electrons. The experimental techniques are so 
refined that one electron at a time can be injected into the 
dot in a fully controllable way. This procedure has shown 
many similarities between quantum dots and atoms, for 
example the existence of shell structure. To emphasize 
these similarities quantum dots are often called artificial 
atoms. The interest in quantum dots is mainly motivated 
by the fact that their properties are tunable through elec- 
trostatic gates and external electric and magnetic fields, 
making these designer atoms promising candidates for 
nanotechnological applications. An additional aspect is 
that quantum dots provide a new type of targets for 
many-body methods. In contrast to atoms they are es- 
sentially two-dimensional and their physical size is sev- 
eral orders of magnitude larger than that of atoms, lead- 
ing e.g. to a much greater sensitivity to magnetic fields. 
Another difference compared to atoms is the strength 
of the overall confinement potential relative to that of 
the electron-electron interaction, which here varies over 
a much wider range. 

The full many-body problem of quantum dots is truly 
complex. A dot is formed when a two-dimensional elec- 
tron gas in an heterostructure layer interface is confined 
also in the xy-plane. The, for this purpose applied gate 
voltage, results in a potential well, the form of which 
is not known. A quantitative account of this trapping 
potential is one of the quantum dot many-body prob- 
lems. Self-consistent solutions of the combined Hartree 
and Poisson equations by Kumar et alX in the early 
nineties indicated that for small particle numbers this 
confining potential is parabolic in shape at least to a 
first approximation. Since then a two-dimensional har- 
monic oscillator potential have been the standard choice 
for studies concentrating on the second many-body prob- 
lem of quantum dots; that of the description of the in- 
teraction among the confined electrons. The efforts to 



give a realistic description of the full physical situation, 
see e.gii~— have, however, underlined that it is important 
to realize the limits of this choice. To start with the 
pure parabolic potential seems to be considerably less 
adequate when the number of electrons is approaching 
twenty. The potential strength is further not indepen- 
dent of the number of electrons put into the dot, an effect 
which is sometimes approximately accounted for by de- 
creasing the strength with the inverse square root of the 
number of electrons^. Finally, the assumption that the 
confining potential is truly two-dimensional is certainly 
an approximation and it will to some extent exaggerate 
the Coulomb repulsion between the electrons. In Ref. 3 
the deviation from the pure two-dimensional situation 
is shown to effectively take the form of an extra poten- 
tial term scaling with the fourth power of the distance 
to the center and which can be both positive and neg- 
ative. Although the deviation is quite small it is found 
that predictions concerning the so called third shell can 
be affected by it. 

There is thus a number of uncertainties in the descrip- 
tion of quantum dots. On the one hand there is the de- 
gree to which real dots deviate from two-dimensionality 
and pure parabolic confinement. On the other hand there 
is the uncertainty in the account of electron correlation 
among the confined electrons. The possible interplay 
among these uncertainties is also an open question. In a 
situation like this it is often an advantage to study one 
problem at a time, since it is then possible to have con- 
trol over the approximations made and quantify their ef- 
fects. We concentrate here on the problem of dot-electron 
correlation. For this we employ a model dot; truly two- 
dimensional, with perfect circular symmetry and with 
a well defined strength of the confining potential. This 
choice is sufficient when the aim is to test the effects of 
the approximations introduced through the approxima- 
tive treatment of the electron-electron interaction. 

Especially the experimental study by Tarucha et. alii 
has worked as a catalyst for a vast number of theoreti- 
cal studies of quantum dots. A review of the theoretical 
efforts until a few years ago has been given by Reimann 
and Manninen 8 . A large number of calculations has been 
done within the framework of Density Functional Theory 
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(DFT)£>£ — and reference therein, but also Hartree-Fock 
(HFJiiriS, Quantum Monte Carlo methods!^ and Con- 
figuration Interaction (CI)^— studies have been per- 
formed. The DFT-studies have been very successful. 
The method obviously accounts for a substantial part 
of the electron-electron interaction. Still, the situation is 
not completely satisfactory since there is no possibility to 
systematically improve the calculations or to estimate the 
size of neglected effects. For just a few electrons the CI- 
approach can produce virtually exact results, provided 
of course that the basis set describes the physical space 
well enough. The size of the full CI problem grows, how- 
ever, very fast with the number of electrons and to the 
best of our knowledge the largest number of electrons in 
a quantum dot studied with CI is six. It would be an 
advantage to also have access to a many-body method 
which introduces only well defined approximations and 
which allows a quantitative estimate of neglected contri- 
butions. The long tradition of accurate calculations in 
atomic physics has shown that Many-Body Perturbation 
Theory (MBPT) has these properties. It is an in princi- 
ple exact method, applicable to any number of electrons, 
and the introduced approximations are precisely defined. 
With MBPT it is possible to start from a good, or even 
reasonable, description of the artificial atom and then re- 
fine this starting point in a controlled way. We are only 
aware of one study on quantum dots that have been done 
with MBPT, the one by Sloggett and Sushko\*i2. They 
did second-order correlation calculations on circular and 
elliptical dots in an environment free of external fields. 

In the present study we use second-order perturbation 
theory to calculate energy spectra for quantum dots with 
and without external magnetic fields. We consider this 
second-order treatment as a first step towards the calcu- 
lation of correlation to high orders through iterative pro- 
cedures, an approach commonly used for atoms^. The 
method is described in Section [TTJ In section lllll we com- 
pare our calculations with experimental results 3,7 , DFT- 
calculations^ and Cl-calculations, our own as well as 
those of Reimann et aL— and discuss the strength and 
limits of the MBPT approach. We have mainly used the 
Hartree-Fock description as starting point for the pertur- 
bation expansion, but we also show examples with a few 
alternative starting points, among them DFT. To obtain 
a complete and finite basis set, well suited to carry out 
the perturbation expansion, we use so called B-splines, 
see e.g. Refi^i. The use of B-splines in atomic physics 
was pioneered by Johnson and Sapirstein22 twenty years 
ago and later it has been the method of choice in a large 
number of studies as reviewed e.g. in Ref<22. We com- 
pare our correlated results to our own HF-calculations, 
thereby highlighting the importance of correlation both 
when the quantum dot is influenced by an external mag- 
netic field and when it is not. We present addition en- 
ergy spectra for the first twelve electrons. The inter- 
esting third shell (electron seven to twelve) is discussed 
in Section HVl Here we investigate several different fill- 
ing sequences and show that correlation effects in many 



cases can change the order of which the shells are filled. 
We note also that the energy of the first excited state 
can be very close to the ground state, in some case less 
than 0.1 meV, which raises the question if it is always the 
ground state which is filled when an additional electron 
is injected in the dot since more than one state may lie 
in the transport window controlled by the source drain 
voltage^. 



II. METHOD 

The essential point in perturbation theory is to divide 
the full Hamiltonian H into a first approximation, h, and 
a correction, U. The first approximation should be easily 
obtainable, in practice it is more or less always chosen to 
be an effective one-particle Hamiltonian, and it should 
describe the system well enough to ensure fast and steady 
convergence of the perturbation expansion. The partition 
is written as 



U. 



(1) 



Here we have chosen to mainly use the Hartree-Fock 
Hamiltonian as h but we have also investigated the pos- 
sibility to use a few other starting points. 

A first approximation to the energy is obtained from 
the expectation value of H , calculated with a wave func- 
tion in the form of a Slater determinant formed from 
eigenstates to h(i). The result is then subsequently re- 
fined through the perturbation expansion. Below we de- 
scribe the calculations step by step. 



A. Single-particle treatment 

For a single particle confined in a circular quantum dot 
the Hamiltonian reads 



h. = 



P 1*22 e 
2m* 2 8m 



-B 2 r 2 + —BL z +g*v b BS z , 
2m* 

(2) 

where B is an external magnetic field applied perpendic- 
ular to the dot. The effective electron mass is denoted 
with m* and the effective g- factor with g* . Throughout 
this work we use m* = 0.067m e and g* = —0.44, corre- 
sponding to bulk values in GaAs. 

The single particle wave functions separate in polar 
coordinates as 

\fnm ems ) = \u nminis {r))\e lm ^)\m s ). (3) 

As discussed in the introduction we expand the radial 
part of our wave functions in so called B-splines labeled 
Bj with coefficients c,-, i.e. 



\u nmema (r)) = }Ci\Bi(r)). 



i=i 



(4) 
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B-splines are piecewise polynomials of a chosen order fc, 
defined on a so called knot sequence and they form a 
complete set in the space defined by the knot sequence 
and the polynomial order—. Here we have typically used 
40 points in the knot sequence, distributed linearly in 
the inner region and then exponentially further out. The 
last knot, defining the box to which we limit our prob- 
lem is around 400 nm. The polynomial order is six 
and combined with the knot sequence this yield 33 ra- 



dial basis functions, 



, (r), for each combination 



(m£,m s ). The lower energy basis functions are physi- 
cal states, while the higher ones are determined mainly 
by the box. The unphysical higher energy states are, 
however, still essential for the completeness of the basis 
set. 

Eqs. ([3]) and (J4| imply that the Schrodinger equation 
can be written as a matrix equation 



where Hi 



He = eBc 



(5) 



(B J e mW \H\B l e mli ) and B. }i 



{B 3 \B^. 



Eq.(EJ) is a generalized eigenvalue problem that can be 
solved with standard numerical routines. The integrals 
in ((HI are calculated with Gaussian quadrature and since 
B-splines are piecewise polynomials this implies that es- 
sentially no numerical error is produced in the integra- 
tion. 



where 



R\ + R\ + (zi - z 2 f 
2R1R2 



(9) 



Assuming a 2D confinement we set Zi = Zo in (|9|). The 
Q m _^ (x)-functions are Legendre functions of the second 
kind and half-integer degree. We evaluate them using a 
modified^ version of software DTORHl.f described in 25 . 

Using (lU) and ([3]) the electron-electron interaction ma- 
trix element ([7]) becomes 



(ab\ — \cd) 
ri2 



4ire r eo 



[u a (ri)u b (rj)\ 



Uc{n)u d (rj)} 



00 
m— — 00 



x(m a s \m c s ){m b s \m d s ll0) 

Note that the angular part of (|T0|) equals zero except if 
m = m a — m c or m = — mj. This is how the degree of 
the Legendre-function in the radial part of (fit)]) is chosen. 
It is also clear from (TTU|) that the electron-electron matrix 
element ([7]) equals zero if state a and c or state b and d 
have different spin directions. 



2. Hartree-Fock 



B. Many-Body treatment 

The next step is to allow for several electrons in the dot 
and then to account for the electron-electron interaction, 
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47re r e | r. t -rj | ' 



(0) 



where e r is the relative dielectric constant which in the 
following calculations is taken to be e r = 12 .4 correspond- 
ing to the bulk value in GaAs. For future convenience we 
define the electron-electron interaction matrix element as 

V 'rV 1 J J 47re r eo|ri - r, | J ' 

(7) 

where a,b,c and d each denote a single quantum state 
i.e. \a) — |n a ,m^,m°). 



1. The Multtpole expansion 

As suggested by Cohl et. al 2 ^, the inverse radial dis- 
tance can be expanded in cylindrical coordinates (i?, <j), z) 
as 
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im(0i-0 2 ) 



(8) 



If the wave function is restricted to be in the form of 
a single Slater determinant, the Hartree-Fock approxi- 
mation can be shown to yield the lowest energy. In this 
approximation each electron is governed by the confining 
potential and an average Hartree-Fock potential formed 
by the other electrons. To account for the latter the 
Hamiltonian matrix H in Eq. ([5]) is modified by the ad- 
dition of a term: 



uf/ = (Bj\uHp\Bi) = 



E 

a<N 



[Bjal^-lBi^-iBjal^-laBi) 
ri2 ri2 

(11) 

The sum here runs over all occupied orbitals, a, defined 
by quantum numbers n, mi, and m s . Eq. ([5J is then 
solved iteratively yielding new and better wave func- 
tions in each iteration until the energies become self- 
consistent. The hereby obtained solution is often labeled 
the unrestricted Hartree-Fock approximation since no ex- 
tra constraints are imposed on u HF . 

One property of the unrestricted Hartrce Fock approx- 
imation deserves special attention. Consider the effects 
of the Hartree-Fock potential on an electron in orbital b, 

(b\u HF \b) = V (6o|^-|6o) - (ba\^-\ab), (12) 

< 4 Ti o Tin 



i<JV 



ri2 



r\2 



where the last term in Eq. (|12l) . the exchange term, is 
non-zero only if orbital a and b have the same spin. For 
configurations where not all electron spins are paired elec- 
trons with the same quantum numbers n and mi, but 
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with different spin directions, will experience different 
potentials. This is in accordance with the physical situ- 
ation, but has also an undesired consequence; the total 
spin, S 2 = s i) 2 > does not commute with the Hartree- 
Fock Hamiltonian. This means that the state vector con- 
structed as a single Slater determinant of Hartree-Fock 
orbitals will not generally be a spin eigenstate. However, 
the full Hamiltonian, Eq. [TJ still commutes with S 2 and 
during the perturbation expansion the spin will eventu- 
ally be restored, provided of course that the perturbation 
expansion converges. Since, in contrast to the energy, the 
total spin of a state is usually known, the expectation 
value of the total spin, (S 2 ), can be used as a measure 
of how converged the perturbation expansion is. It can 
also be used as an indication of when the Hartree-Fock 
description is too far away from the physical situation 
to be a good enough starting point. This is discussed 
further in Sections Mil and IIVI 



3. second-order correction to a Hartree-Fock starting point 

The leading energy correction to the Hartree-Fock 
starting point is of second order in the perturbation 
(defined in Eq. ((TJ). When h = h s + uhf and U — 
J2i<j ~F~ ~ the corresponding corrections 

to the wave function will not include any single excita- 
tions. This is usually referred to as Brillouin's theorem 
and is discussed in standard Many-Body theory text- 
books, see e.g. Lindgren and Morrison^. Starting from 
the HF-Hamiltonian for N electrons in the dot we write 
the second-order correction to the energy 



SE 



(2) 



x - v \(rs\±\ab)\ 2 -(ba\±\rs)(r S \±\ab) 

a<b<Nr,s>N ^ + ^ ^ 

r^s 

(13) 

where thus both r and s are unoccupied states. 

Since B-splines are used for the expansion of the radial 
part of the wave functions there is a finite number of 
radial quantum numbers (n) to sum over in the second 
sum of Eq. (|13[) . However, in principle there is still an 
infinite number of angular quantum numbers (mi) to sum 
over in the same sum. In praxis this summation has to 
be truncated and the effects of this truncation will be 
discussed in section! 



Hartree-Fock) an alternative starting point might con- 
verge faster, or even provide convergence in regions where 
it cannot be achieved with Hartree-Fock. First of all we 
start with the simplest possible starting point; the pure 
one-electron wave functions. In this case the basis set 
consists of the solutions to the pure 2D harmonic oscil- 
lator in the chosen box and we treat the whole electron- 
electron interaction as the perturbation. The second al- 
ternative starting point is the Local Density Approxima- 
tion (LDA) . That is we switch the second term in Eq. (jTTj) 

to a(Bj\4a* B ^J 2p ^ \Bi) , where p{r) is the radial electron 
density and a is called Slaters exchange parameter and 
is usually set to one. Both these starting Hamiltonians 
are defined with only local potentials and will thus com- 
mute with the total spin. The third alternative starting 
point is a reduced exchange HF, i.e. the exchange term 
(the second term) in Eq. (|11[) is simply multiplied with 
a constant a < 1. When using these alternative start- 
ing points, one must in contrast to the Hartree-Fock case 
include single excitations in the perturbation expansion. 

The second-order perturbation correction then be- 
comes 



SE 



V= E E 

a,b<N r>N 



\{r\V ex \a) - (rb\±\ba)f 



E E 

a<b<N r,s>N 
r^s 



\( rs \±\ab)\ 2 - (ba\±\rs)(r S \£- 2 \ab) 



(14) 



where V ex is the chosen exchange operator. From this 
expression it is also clear that the first term yields zero 
in the pure Hartree-Fock case, i.e. then all single excita- 
tions cancel. 



5. Full CI treatment of the two body problem 

To investigate how well second-order many-body per- 
turbation theory performs we have for the simple case of 
two electrons also solved the full CI problem. We then 
diagonalize the matrix that consists of all the elements 
of the form 



1 



Hji = {mn\. h\ + h 2 s + — \op) 



ri2 



(15) 



4- Other starting points than Hartree-Fock 



for given values of Ml = ^ mi and Ms = m s of 
our electron pairs {|mn)i}. Following the selection rules 
produced by Eq. (|10[) we get the conditions m£ + = 
mT + mj, m™ = m° and m™ = m?. 



In principle any starting point, with wave functions 
close enough to the true wave functions (to ensure con- 
vergence of the perturbation expansion), can work as a 
starting point for MBPT. We have in addition to HF in- 
vestigated three alternative starting points. If there are 
important cancellations between the full exchange (in- 
cluded in Hartree-Fock) and correlation (not included in 



III. VALIDATION OF THE METHOD 

The main purpose of this work is to investigate the 
usability of many body perturbation theory on (GaAs) 
quantum dots. Therfore we have in this section compared 
our results with results from other theoretical works. 
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FIG. 1: Second-order perturbation theory correction to the 
energy as function of max(n) (squares) and max(|m^|) (cir- 
cles) in the second sum of Eq. (I13p for the two electron dot 
with the confinement strength ftw = 6 meV. Note that the 
sum over mi converges faster than the sum over n. 



Our energies are generally given in meV. For easy 
comparison with other calculations it should be noted 
that the scaled atomic unit for energy is 1 Hartree* = 
1 Hartree (m*/(m e e^)) s» 11.857 meV, with to* = 
0.067m e and e r = 12.4. The scaled Bohr radius is 
a* B = (e r m e /m*) clb ~ 9.794 nm. 



A. The two electron dot 

Fig. [T] shows the second-order many-body perturba- 
tion correction to the energy, Eq. (IT3|) . as function of 
max(n) (squares) and max(|m^|) (circles) respectively for 
the two electron dot with ftui — 6 meV. It clearly illus- 
trates that both curves converge but also that the sum 
over uig converges faster than the sum over n. Due to 
this we have throughout our calculations used all radial 
basis functions and as many angular basis functions that 
are needed for convergence. One should, however, notice 
that the relative convergence as a function of max(n) and 
max(|m^|) varies with the confinement strength and oc- 
cupation number. Weak potentials (ftw < 2 meV) usually 
produce the opposite picture i.e. a faster convergence for 
n than for mg. For confinement strengths (tvjj > 3 meV) 
and most occupation numbers the trend shown in Fig. ([1]) 
is, however, typical. 



1. Comparison between different starting points 

In Fig. [2] a) and b) comparisons between HF (i.e. the 
expectation value of the full Hamiltonian with a Slater 
determinant of Hartree-Fock orbitals, labelled HF + 1st 
order MBPT in Figs. [5] and [3]), the alternative starting 
points discussed in section III B 41 second-order MBPT 
(HF or alternative starting point + second-order correla- 
tion) and CI calculations for the ground state in the two 
electron dot are made. Both the second-order MBPT 
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FIG. 2: The quotient between the calculated energies (of re- 
spective method) and the Cl-energy as functions of the con- 
finement strength, hui, for the ground state in the two electron 
dot. In a) the results from HF, a reduced exchange HF with 
a = 0.7 and 2nd order MBPT using wave functions from 
respective method are plotted. In b) the results from LDA- 
calculations (with two different alphas) + 1st and 2nd order 
MBPT are plotted. For reference the results from calcula- 
tions where we have used the one-electron wave functions as 
a starting point for the perturbation expansion (taking the 
whole electron-electron interaction as the perturbation) have 
been plotted in both a) and b). 



and CI results have been produced with all radial ba- 
sis functions (33 for each combination of me and m s ) 
and — 6 < mi < 6. It is clear from Fig. [5^) that second- 
order correlation here is the dominating correction to the 
Hartree-Fock result. Even for ftu = 2 meV the differ- 
ence compared to the CI result decreases with one order 
of magnitude when it is included. For stronger confine- 
ments the difference to CI is hardly visible. As expected, 
the performance of both HF and second-order MBPT is 
improved when stronger confinement strengths are con- 
sidered. For the weakest confinement strength calculated 
here (ftio — 1 meV) the pure Hartree-Fock approxima- 
tion gives unphysical wave functions in the sense that 
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the spin up and the spin down wave functions differ, re- 
sulting in a non-zero (S 2 ). This shows up in figure a) as 
a broken trend (all of a sudden an overestimation of the 
energy instead of an underestimation) for the pure HF 
+ second-order correlation curve at Hu) = 1 meV. For 
all other potential strengths (S 2 ) is zero to well below 
the numerical precision 10 -6 ) for both the Hartree- 
Fock and second-order MBPT wave functions, while for 
the fvjj = 1 meV calculation (S* 2 ) = 0.33 and 0.26 for 
the Hartree-Fock and second-order MBPT calculations 
respectively. It should be noted that at tkJ = 1 meV the 
energy of the second-order MBPT calculation is still only 
4% larger than the CFvalue (although the wave functions 
are unphysical) and that probably the state will converge 
to (S 2 ) = when MBPT is performed to all orders. All 
other tested starting points yield (S 2 ) = for this con- 
finement strength, but still their energy estimates after 
second-order MBPT are worse. This shows that con- 
served spin does not necessary yield good energies and 
broken spin symmetry does not necessary yield bad en- 
ergy estimations. We note that the reduced exchange 
Hartree-Fock, displayed in Fig. [5^) , seems to be a fruit- 
ful starting point for perturbation theory although the 
results after second-order are slightly worse than after 
the full exchange Hartree-Fock + second-order MBPT. 
For Huj = 1 meV the reduced exchange HF with a = 0.7 
still gives (S 2 ) — 0, i.e. the onset of spin contamination 
is delayed on the expense of proximity to the Cl-energy. 
To put it in another way, if we lower a, the corresponding 
curve in Fig [5] a) will be lower (and thus further from the 
correct Cl-curve) , but the spin contamination onset will 
appear for a weaker confinement strength. This freedom 
could be useful when doing MBPT to all orders. 

From Fig. [2] b) we conclude that LDA with a > 1 
might be a useful starting point for perturbation theory 
calculations to all orders but not a good option for 2nd 
order calculations, at least not for weak potentials. LDA 
calculations with a = 1, however, seems to be a bad 
starting point for MBPT, at least for the two electron 
case, since it gives almost identical results after second 
order as using the pure one electron wave functions as 
starting point. LDA might still work better as a starting 
point when more electrons are added to the dot. 

Finally the comparison with the pure one-electron 
wave functions in Fig. [2] clearly illustrates how much of 
an improvement it is to start the perturbation expan- 
sion from wave functions that already include some of 
the electron-electron interaction, especially for weaker 
potentials. This becomes even more clear in Fig.[3]where 
we present the results from Fig.[5]a) in another way. Here 
we have plotted Euethod I Eq as functions of Iq / a* B where 
Eq = Huj is the single particle energy and Iq — ^Jh/(m*L>S) 
is the characteristic length of the dot. It demonstrates 
what an extraordinary improvement it is to start from 
Hartree-Fock compared to starting with the one-electron 
wave functions when doing second-order MBPT for low 
electron densities (high lo/a* B ). It also seems as there 
is a region where the Hartree-Fock starting point would 



-A- no e-e interaction 
-A-no e-e basis + 1st order MBPT 
no e-e basis + 2nd order MBPT 
-V-HF°"°- 7 + Istord MBPT 
-V-HF"" ' 7 + 2nd order MBPT 
-0-HF + Istord MBPT 
-0-HF + 2nd order MBPT 




FIG. 3: E/Eo for the two-electron dot and for different meth- 
ods as functions of lo /a B , where lo = yhj (m*ui) is the char- 
acteristic dot length, Eo = ftw is the single particle energy, 
and a* B is the effective Bohr radius. Small values of lo/a* B 
correspond to stronger confinement and therefore a faster ex- 
pected convergence rate of a perturbation expansion. 



yield a convergent perturbation expansion while taking 
the whole electron-electron interaction as the perturba- 
tion would not. 



B. The six electron dot 

In Fig. [4l a comparison between our HF and second- 
order MBPT calculations on the ground state of the six 
electron dot is made with a DFT calculation in the Local 
Spin Density Approximation(LSDA) as well as with a CI 
calculation, both by Reimann et al*£. They performed 
their calculations for seven different electron densities 
here translated to potential strengths. Let us first focus 
on the results for the two highest densities, correspond- 
ing to a Wigner-Seitz radius r s = l.Oa^ and r s = 1.5a B 
which translates to confinement strengths of fiw w 7.58 
meV and rs 4.12 meV respectively. The reason that we 
want to separate the comparison for those confinement 
strengths is that our Hartree-Fock calculations yield solu- 
tions with (S 2 ) > for the weaker confinement strengths. 
A similar behavior was seen by Sloggett et al^ in their 
unrestricted HF calculations. Therefore the results for 
the weaker potentials overestimate the energy in a un- 
physical manner; compare the above discussion around 
Fig. [2] a). The Cl-method however always yield (S 2 ) = 
for the closed 6 electron shell and consequently a compar- 
ison with spin contaminated results would here, in some 
sense, be misleading. It should be emphasized that the 
spin contamination is a feature of our choice of starting 
point and not a problem with MBPT in itself. 

To make comparison easy all energies are normalized 
to the corresponding Cl-value. The figure clearly illus- 
trates, for the two stronger confinement strengths, that 
while the HF results overshoot the CI energy by between 
3.5% and 4.5% the second order MBPT calculations im- 
prove the results significantly. Already for max(|m^|) = 1 
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FIG. 4: Comparison between our HF and second-order 
MBPT results for the six electron dot in the ground state with 
Ml OT = and ST OT = 0, with the LSDA and CI calculations 
by Reimann et al 1(i . The second-order MBPT calculations 
include the full sum over the complete radial basis set (corre- 
sponding to all n-values) and with max(|m£|) = 1, 2, 3, . . . , 30 
for the two strongest potentials. For clarity only the curves 
with max(|m^|) = 1,4,6 and 30 have been labeled. The HF 
and the second-order MBPT with max(|m^|) = 30 curves are 
plotted for all potential strengths calculated by Reimann et 
al. Moreover, the values of {S 2 } for the HF and the second- 
order MBPT with max(|m^|) = 30 have been plotted in the 
figure. 



the energy only overshoots the CI value with between 
2.5% and 3.5% while the second-order MBPT energy for 
max(m£) = 4 is almost spot on the CI energy. However, 
with max(m<) = 30 the second-order MBPT gives some- 
where between 0.5% and 1% lower energy than the CI 
calculation. We note that the CI calculation by Reimann 
et al. was made with a truncated basis set consisting of 
the states occupying the eight lowest harmonic oscilla- 
tor shells. This means e.g. that their basis set includes 
only two states with (\mi |) = 5 and one with = 6. 

Within this space all possible six electron determinants 
were formed. After neglecting some determinants with a 
total energy larger than a chosen cutoff, the Hamiltonian 
matrix was constructed and diagonalized. Fig. 2] indi- 
cates that the basis set used in Ref J£ was not saturated 
to the extent probed here, since almost all interactions 
with \rri(\ > 4 were neglected. According to Reimann 
et al. they used a maximum of 108 375 Slater deter- 
minants while we, through perturbation theory, use a 



maximum of 980 366 Slater determinants. The differ- 
ence of our max(|m^|) = 30 results and their CI results 
are thus not unreasonable. Since Reimann et al. solved 
the full CI problem, the matrix to diagonalize is huge 
and it is, according to the authors, not feasible to use 
an even larger basis set. An alternative could be to in- 
clude more basis functions, but restrict the excitations 
to single, doubles and perhaps triples. The domination 
of double excitations is well established in atomic cal- 
culations, see e.g. the discussion in Refill. It should 
however be noted that the difference between the results 
concerns the fine details. Our converged results are less 
than one percent lower than those of Reimann et al. and 
when using approximately the same basis set as they did 
(max(|m^ |) = 4) the difference between the results is vir- 
tually zero. Moreover, we see for the two strongest po- 
tentials the same trend as we saw in the two-electron 
case, namely that the HF, MBPT and CI results tend 
towards one another with increasing potential strength. 
This trend is not seen for the LSDA approach. 

Finally, Fig. [4] shows, for the five weakest poten- 
tials, that our HF results get increasingly spin contami- 
nated when the potential is weakened. Hereby the HF- 
approximation artificially lowers its energy and subse- 
quently this leads to an overestimation of the second- 
order MBPT energies for these potential strengths. Sur- 
prisingly, however, the energy is never more than just 
above 2% over the Cl-results even when (S 2 ) > 2. Note 
also that MBPT improves the HF-value of (S 2 ) as it 
should. 



C. Correlation in an external Magnetic Field 

The behavior of quantum dots in an external magnetic 
field applied perpendicular to the dot has previously been 
examined many times both experimentally e.gi^ ' 24 ' 28 and 
theoretically e.g£ ' 18 i 29 . The chemical potentials /J.(N) = 
E(N) — E(N—1) plotted versus the magnetic field usually 
show a rich structure, including e.g. state switching and 
occupation of the lowest Landau band at high magnetic 
fields. 

Fig J5] shows the chemical potentials for N = 1,2, ... ,6 
as functions of the magnetic field according to our HF 
(dashed curves) and second-order MBPT with — 10 < 
mi < 10 (full curves) calculations for the potential 
strength tiui ~ 5 meV. We have here limited ourselves 
to the first six chemical potentials calculated at selected 
magnetic field strengths (shown by the marks in the fig- 
ure). We emphasize again that our intention here is 
rather to test the capability of MBPT in the field of quan- 
tum dots than to provide a true description of the whole 
experimental situation. With increasing particle number 
MBPT naturally becomes more cumbersome, but mag- 
netic field calculations are feasible at least up to N = 20. 

First note the significant difference between the HF 
and second-order MBPT results. Once again correlation 
proves to be extremely important in circular quantum 
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FIG. 5: The chemical potentials for TV = 1,2,. ..,6 as func- 
tions of the external magnetic field according to HF (dashed 
curve) and second-order MBPT (full curve) calculations for 
the potential strength hu — 5 meV. Note the big difference 
between the two different models regarding the behavior of 
the chemical potentials when the magnetic field varies. 



FIG. 6: The ground state addition energy spectra for dots 
with Hlj = 5 meV (a) and hw — 7 meV (b). The squares 
(circles) represent the addition energy spectra according to 
HF (second-order MBPT). It is clear that the second-order 
MBPT-spectra imply closer resemblances to the experimental 
picture in Tarucha? than the HF-spectrum. 



IV. RESULTS 



dots. With our correlated results we also note a close re- 
semblance both to the experimental work by Tarucha et 
alZ and to the current spin-density calculation by Stef- 
fens et al*&, made with the same potential and material 
parameters as used here. (Note that Ref^ defines the 
chemical potentials as fi(N) = E(N + 1) - E(N), shift- 
ing all curves one unit in N). An example of the impor- 
tance of correlation is the four-electron dot that switches 
state from | m, \M L \, S) = |0, 0, 1) to |0, 2, 0) at ap- 
proximately 1 T in the HF calculations and at approx- 
imately 0.2 T in the correlated calculations. We want 
to emphasize that we have found the exact position of 
this switch to be very sensitive to the potential strength 
and to the value of g*. The big difference concerning 
the magnetic field where this switch occurs can proba- 
bly be attributed to the HF tendency to strongly favor 
spin- alignment. This is an effect originating from the in- 
clusion of full exchange, but no correlation. Inclusion of 
second-order correlation energy cures this problem. Fi- 
nally we note that the N — 5 switch from |0, 1, |) to 
|0,4, i) in our correlated calculations takes place some- 
where around 1.2 T which is also in agreement with both 
mentioned studies. 



A. The addition energy spectra 

The so called addition energy spectra, with the addi- 
tion energy denned as A(N) = E(N + 1) - 2E(N) + 
E(N — 1), have been widely used to illustrate the shell 
structure in quantum dots. Main peaks at N = 2, 6, 12 
and 20, indicating closed shells, and subpeaks at N = 4, 9 
and 16, due to maximized spin at half filled shells, have 
been interpreted as the signature for truly circular quan- 
tum dots^. Experimental deviations from this behavior 
have been interpreted as being due to nonparabolicitics 
of the confining potential or due to 3D-effects 3 -. We here 
show that correlation effects in a true 2D harmonic po- 
tential can in fact generate an addition energy spectrum 
with similar deviations. 

In this work we limit ourselves to the first three shells 
since it seems as the experimental situation is such that 
the validity of the 2D harmonic oscillator model becomes 
questionable with increasing particle number 3 . Calcula- 
tions of dots with larger N could, however, readily be 
made with our procedure. The addition energy spectra 
are produced with —10 < mi < 10. The filling order 
for the first six electrons is straight forward. When the 
seventh electron is added to the dot the third shell starts 



9 



to fill. With a pure circular harmonic oscillator potential 
and no electron-electron interaction the | 0, ±2, ±1) and 
| 1,0, ±5) one-particle states are completely degenerate. 
This degeneracy is lifted by the electron-electron inter- 
action, but not more than that the energies have to be 
studied in detail in order to determine the filling order. 
Similar conclusions, that the filling order is very sensitive 
to small perturbations, have been drawn by Matagne et 
al&, who studied the influence of non-harmonic 3D ef- 
fects. Our focus is instead the detailed description of the 
electron-electron interaction. For N = 7 — 11 we have 
thus calculated all third shell configurations, and for each 
configuration considered the maximum spin. The results 
are found in Tablefl] For each number of electrons we can 
identify a ground state, which sometimes differs between 
HF and MBPT. These ground states are used when cre- 
ating Fig. [n^) and b). The energy gap to the first excited 
state is sometimes very small and the possibility of alter- 
native filling orders will be discussed in the next section. 

Fig. [Bh) and Fig. [5b) thus show the ground state ad- 
dition energy spectra up to N — 12 according to the 
Hartree-Fock model as well as to second-order MBPT 
for haj — 5 and 7 meV. Note first the big difference be- 
tween the HF and MBPT spectra. These figures clearly 
illustrate how important correlation effects are in these 
systems. Admittedly the HF-spectra show peaks at 
N = 4,6,9 and 12 but the relative size of the addition 
energy between closed and half-filled shells is not con- 
sistent with the experimental picture 3 '-. The second- 
order MBPT-spectra have in contrast clear main peaks 
at N = 2, 6 and 12, indicating closed shells, and a A = 4 
subpeak indicating maximized spin for the half filled 
shell. For the fku — 7 meV spectrum the subpeak at 
N — 9 is also clear but for the hiu = 5 meV spectrum the 
subpeak at N — 9 is substituted by subpeaks at N = 8 
and 10. The behavior of the addition energy spectra in 
this, the third shell, will be discussed in detail below. 



1. Filling of the third shell 

The filling of the third shell has previously been exam- 
ined by Matagne et al£ both experimentally and theo- 
retically. In their theoretical description they use a 3D 
DFT model with the possibility to introduce a nonhar- 
monic perturbation that can change the ground states 
in the third shell and thereby alter the addition energy 
spectra. They then compare their theoretical description 
with different experimental addition energy spectra and 
argue how large deviation from the circular shape they 
have in the different experimental setups. They conclude 
that a clear dip at TV = 7 followed by a peak at N = 8 
or 9 is a signature of maximized spin at half filled shell 



and that a dip at N = 7 and the filling sequence 

N N \ 

5>M£m||,s)=|0,2,-)^| 0,0,1)^ 

i=7 i=7 / 

I 1,0, |) =>| 1,2,1) =H 1,0, i) =H 2,0,0) (16) 

for the six electrons to enter the third shell is a signature 
of a "near ideal artificial atom" . This is also the filling 
sequence we find using the HF- approximation. As seen 
in Fig. [5^) and b) there is then indeed also a dip at 
N = 7 and a peak at N = 9. The dip at N = 7 is 
further supported by the DFT calculation by Reimann 
et al£Q. In contrast the experiment by Tarucha et alZ 
did not show the N — 7 dip. In Ref4 this is explained 
by deviations from circular symmetry for the specific dot 
used in Ref* 7 -. As will be seen below our many-body 
calculations give in several cases different ground states 
and thus favor a different filling order than Eq|T6l 

TableUshows the ground state and excited states ener- 
gies of the third shell according to HF and second-order 
MBPT for hw = 5 meV and frw = 7 meV. Notice that 
the different methods yield different ground states for 
the 8, 10 and 11-electron systems although both poten- 
tial strengths yield the same ground states. Note also the 
small excitation gap between the correlated ground and 
first excited state that occurs in some cases. For example 
between the | 0,2,4) and | 1,0, 4) seven-electron states 
in the huj = 5 meV dot the energy difference is 0.07 meV, 
between the | 0,0, 1) and | 0,4,0) eight-electron state in 
the huj = 7 meV dot the energy difference is 0.12 meV 
and between the | 1,0,^) and | 2,2,^) eleven-electron 
states in the hu> — 7 meV dot the energy difference is 
only 0.04 meV. The (1, 0, §) state at N = 9 seems, how- 
ever, relatively stable for both potential strengths with 
excitation gaps of 0.41 and 0.54 meV. Surprisingly for 
both the hw — 5 and 7 meV the calculations including 
correlations indicate the ground state third shell filling 
sequence 

I 0,2,i)=H 0,4,0) =H 1,0,|> => 

I 0,0,0) H 2,2,^) =H 2,0,0) (17) 

for N = 7 — 12. Note that this sequence implies a spin- 
flip of the electrons already in the dot when the ninth 
and tenth electrons are added. Only the seven-electron 
dot and the nine-electron dot here have the same ground 
state as in HF (whose filling sequence coincides with that 
preferred in RefS). Matagne et al. also discuss that the 
behavior of the dot examined in Ref. 7 for small magnetic 
fields implies the sequence 

I o,2,i)=H 0,4,0} =>| i,2,i}=* 

I 0,0,0) =H i,o,i)=H 2,0,0), (is) 

but tend to attribute this to deviations from circular 
shape. This filling sequence is indeed much closer to the 
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TABLE I: Energy of the ground and third shell excited state for 7-11 electron dots with ftu = 5 and 7 meV. The notation 
(Y2iLi n > I -Mill "SO t° label the state is used. The ground state energy according to Hartree-Fock (HF energy) and to HF + 
second-order MBPT (Correlated energy) and for respective N and potential strength is marked in bold. 



# e- 






hw = 


5 meV 






hw = 


7 meV 




7 


State 


(0,2, 


k) 


(1,0, 


4) 


(0,2, 


1) 


(1,0, 


*) 




HF energy [meV] 




168.02 


168.67 


215.80 


216.58 




Correlated energy 


[meV] 


162.08 


162.15 


208.96 


209.52 


8 


State 


(0,0,1) 


(0,4, 0) 


(1,2,1) 


(2,0,0) 


(0,0,1) 


(0,4,0) 


(1,2,1) 


(2,0,0) 




HF energy [meV] 




210.69 


212.33 


211.66 


214.00 


270.66 


272.20 


271.51 


274.32 




Correlated energy 


[meV] 


205.23 


204.40 


204.66 


205.02 


263.82 


263.70 


263.85 


264.65 


9 


State 


(i,o,|) 


(0,2,^) 


(1,4, i) 


(2,2,±) 


(1,0, f) 


(0,2,i) 


(1,4, *) 


(2,2,i) 




HF energy [meV] 




257.69 


259.24 


259.28 


260.56 


330.25 


332.17 


332.14 


333.64 




Correlated energy 


[meV] 


250.54 


251.35 


250.95 


251.00 


322.27 


322.81 


323.06 


323.37 


10 


State 


(1,2,1) 


(0,0,0) 


(2,0,1) 


(2,4,0) 


(1,2,1) 


(0,0,0) 


(2,0,1) 


(2,4,0) 




HF energy [meV] 




309.27 


310.64 


310.17 


311.06 


395.72 


397.20 


396.73 


397.78 




Correlated energy 


[meV] 


300.49 


300.00 


300.25 


300.52 


385.92 


385.76 


386.06 


386.49 


11 


State 


(1,0, 


*) 


(2,2, 


i) 


(1,0, 


4) 


(2,2, 






HF energy [meV] 




363.72 


364.49 


464.77 


465.57 




Correlated energy 


[meV] 


353.66 


353.19 


453.47 


453.43 



TABLE II: Expectation values of S 2 for the cases where correlation switches ground states in the third shell. The state labeled 
"Ground State" is the ground state according to second-order MBPT while the state labeled "Excited State" is the ground 
state according to Hartree-Fock but an excited state according to second-order MBPT. 







hw 


= 5meV 






hjj - 


= 7meV 






Ground State 


Excited state 


Ground State 


Excited state 




E [meV] 


(S 2 ) 


E [meV] 


(S 2 ) 


E [meV] 




E [meV] 


(S 2 ) 


8 


HF 


212.33 


0.00 


210.69 


2.70 


272.20 


0.00 


270.66 


2.30 




2 nd -ord MBPT 


204.40 


0.00 


205.23 


2.58 


263.70 


0.00 


263.82 


2.22 




Exact 









2 









2 


10 


HF 


310.64 


0.00 


309.27 


2.21 


397.20 


0.00 


395.72 


2.08 




2 nd -ord MBPT 


300.00 


0.00 


300.49 


2.15 


385.76 


0.00 


385.92 


2.05 




Exact 
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11 


HF 


364.49 


0.77 


363.72 


0.99 


465.57 


0.758 


464.77 


0.82 




2 nd -ord MBPT 


353.19 


0.76 


353.66 


0.93 


453.43 


0.755 


453.47 


0.79 




Exact 




0.75 




0.75 




0.75 




0.75 



ground states we have obtained with a perfect circular 
potential. This indicates the possibility that many-body 
effects usually neglected could have an effect similar to 
that of imperfections in the dot construction. We note 
in passing that Sloggett and Sushko\*i2 support our find- 
ing of a spin-zero ground-state for ten electrons, although 
their calculation was done with a stronger potential. The 
different configurations for nine electrons in Eq. [T7] and 
Eq. [18] can be due to the fact that the experimental situ- 
ation favors population of an excited state since popula- 
tion of the ground state would require a spin flip. How- 
ever, if we produce a spectrum with this filling sequence, 
we get a large dip at N = 9. Similarly, when the eleventh 
electron is injected, the population of our ground state 
would require a configuration change of the electrons al- 
ready in the dot. 

In Fig. [7] and Fig. [5] addition energy spectra are shown 
assuming different filling orders for 5 meV and 7 meV, 
respectively. In each figure the calculated ground state 



filling sequence is shown in the uppermost panel, labeled 
a), and then the other panels, e) - f), show selected ex- 
cited state filling sequences. Note that even though the 
same filling sequences are used in Fig. [7J and Fig. [5J the 
addition energy spectra differ between these rather close 
potential strengths. We can thus conclude that a given 
filling sequence does not yield a unique addition energy 
spectra since the relative energies of the ground and ex- 
cited states are very sensitive to the exact form of the 
potential. Furthermore we agree with Matagne et al& 
that full spin alignment for the nine-electron ground state 
does not guarantee a peak in the addition energy spec- 
trum as seen in Fig. [Jja) and b). Moreover we see that 
the spectra that resemble the experimental one in Fig. 
3a) of Refi^ (a clear dip at N — 7 and 10 and a clear 
peak at N = 9) are Fig. [Tb) and Fig. [5b). Finally we 
see that Fig. [Sfc) resembles the experimental situation in 
RefZ- (dips at N = 8 and 10 with a peak at N = 9) 
the most. We certainly do not claim that these filling 
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FIG. 7: Ground state, a), and selected excited state, b)-e), 
addition energy spectra for ftw = 5 meV according to second- 
order MBPT. The notation (£? =7 n\ | E l = 7 m\\,S) to label 
the states is used. Note the big differences between the differ- 
ent spectra. For example the ground state spectrum, a), has 
peaks at TV = 8, 10 while spectrum b) has a peak at TV = 8 
and the rest have a peak at TV — 9. That is, even if the spin 
is maximized at half filled shell (TV = 9) there is not always 
a peak there as seen in subfigure a) and b). Subfigure e) re- 
sembles the experimental results of Ref.— best with dips at 
TV — 7 and 10 and a peak at TV — 9. Moreover, combining the 
addition energies for TV = 6, 7, 8 of sequence c) or d) with the 
addition energies for TV = 10, 11, 12 of sequence e) would give 
a spectrum that closely resembles the experimental situation 
in Refj£ with dips at TV = 8 and 10 and a peak at TV = 9. 



FIG. 8: Ground state, a), and selected excited state, b)-e), 
addition energy spectra for ftw — 7 meV according to second- 
order MBPT. The notation n% | J2i=r m H S ) to label 
the states is used. Note the big differences between the differ- 
ent spectra. Note also that all spectra have peaks at TV = 9. 
Even though the ground state spectra for N = 7,8 and 9 re- 
semble the experimental results of Ref.—, the dip at TV = 11 is 
uncharacteristic when compared with the experimental results 
of Ref. 7 and Ref. 3 . Subfi gure b) resembles the experimental 
result in Ref.— the most with a peak at TV = 9 and dips at 
TV = 7 and 10 while subfigure c) resembles the experimental 
results of Ref.— the most with a peak at TV = 9 and dips at 
TV = 8 and 10. 
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FIG. 9: {S 2 } according to Hartree-Fock and second-order 
MBPT calculations as functions of the potential strength for 
the 7 electron ground and excited state. 



overestimated with more than a couple of percent due to 
spin contamination. For the excited state the spin con- 
tamination is so small (for the 5 and 7 meV calculations) 
that it should not affect the conclusions from this work. 
Moreover we see that, as expected, correlation improves 
the value of (S 2 ). 

Table [TT] presents the spin contamination for the sys- 
tems in the third shell where correlation switched the 
ground state, namely the 8, 10 and 11 electron systems. 
We see that the ground states, according to our corre- 
lated results, are not spin contaminated to any relevant 
magnitude. All the excited states are however spin con- 
taminated. As shown in Fig. 2J spin contamination can 
lower the HF energy and raise the second-order MBPT 
energy. The ground state energy switches could thus be 
an artifact of our starting point. Energy wise however 
the correlated energies should lie much closer to the true 
values than the HF-energies. 



V. CONCLUSIONS 



sequences are those really obtained in the mentioned ex- 
periments. However, we want to stress that great care 
must be taken when conclusions are drawn from com- 
parisons between theoretical and experimental addition 
energy spectra. 

2. Spin contamination in the third shell 

Fig. IH1 shows the expectation value of the total spin, 
(S 2 ), according to Hartree-Fock and second-order MBPT 
calculations as functions of the potential strength for the 
7 electron ground and excited state. The figure depicts 
the drastic onset of spin contamination for weak poten- 
tials. While especially the correlated results, but also 
the HF-results, converge towards the correct value for 
potentials > 10 meV the situation is worse for weaker 
potentials. We see that for the ground state the exam- 
ined confinement strengths in this article (huj — 5 or 7 
meV) lie on the onset of the spin density wave. It is 
hard to say how much this spin contamination affects 
the energy values but when compared with the conclu- 
sions drawn from Fig. [2] and 01 the energy should not be 



We have shown that the addition of second-order cor- 
relation improves the Hartree-Fock description of two- 
dimensional few-electron quantum dots significantly. Our 
results indicate that details in the addition energy spec- 
tra often attributed to 3D-effects or deviations from cir- 
cular symmetry, are indeed sensitive to the detailed de- 
scription of electron correlation on more or less the same 
level. Without precise knowledge of the many-body ef- 
fects far reaching conclusions about dot properties from 
the addition energy spectra might not be correct. 

As a next step we want to include pair-correlation to 
higher orders to be able to determine energies with quan- 
titative errors below O.lmeV. We will then use several 
different starting potentials to be able to address also 
weak confining potentials where the Hartree-Fock start- 
ing point fails. 
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